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ABSTRACT 

We present transmission spectra of the hot Jupiter HD 189733b taken with the 



Space Telescope Imaging Spectrograph aboard HST. The spectra cover the wavelength 
range 5808 - 6380 A with a resolving power of R — 5000. We detect absorption 
from the Na I doublet within the exoplanet's atmosphere at the 9a confidence level 
within a 5 A band (absorption depth 0.09 ± 0.01 %) and use the data to measure the 
doublet's spectral absorption profile. We detect only the narrow cores of the doublet. 
The narrowness of the feature could be due to an obscuring high-altitude haze of an 
unknown composition or a significantly sub-solar Na I abundance hiding the line wings 
beneath a H_2 Rayleigh signature. These observations are consistent with previous 
broad-band spectroscopy from ACS and STIS, where a featureless spectrum was seen. 
CN | We also investigate the effects of starspots on the Na I line profile, finding that their 

impact is minimal and within errors in the sodium feature. 

We compare the spectral absorption profile over 5.5 scale heights with model spec- 
tral absorption profiles and constrain the temperature at different atmospheric regions, 
allowing us to construct a vertical temperature profile. We identify two temperature 
regimes; a 1280 ± 240 K region derived from the Na I doublet line wings correspond- 
ing to altitudes below ~ 500 km, and a 2800 ± 400 K region derived from the Na I 
^ ■ doublet line cores corresponding to altitudes from ~ 500 — 4000 km. The zero altitude 

is defined by the white-light radius of Rp/R± = 0.15628 ± 0.00009. The temperature 
rises with altitude, which is likely evidence of a thermosphere. 

The absolute pressure scale depends on the species responsible for the Rayleigh 
signature and its abundance. We discuss a plausible scenario for this species, a high- 
altitude silicate haze, and the atmospheric temperature-pressure profile that results. 
In this case, the high altitude temperature rise for HD 189733b occurs at pressures of 
10-5-lCT 8 bar. 

Key words: techniques: spectroscopic - planetary systems - stars: individual: 
HD 189733 



1 INTRODUCTION 

Transiting planets allow us to study exoplanet atmospheres 
through transmission spectroscopy, as the light from the host 
star passes through the planet's atmosphere. Hot Jupiters 
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offer the best opportunity for these studies due to their large 
size and often large atmospheric scale heights. 

Theory predicts the p resence of strong Na I D lines 
in ho t Jupiter atmospheres (|Seager fc S assclov 2 0001 ; iBrownl 
l200ll) . and observations have so far agreed, with atomic 
sodium detected in hot Jupiters dCharbonneau et al.l 12002 ; 
iRedfield et all 120081 : ISnellen et al.ll2008l : ISing et alj |2008a: 
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I Wood et al.ll201ll ; |j ensen et all 1201 lh . However, the spec- 
tral shape of the sodium doublet absorption has been mea- 
sured only for very few planets. The doublet samples a 
large range of altitudes in the atmosphere and by measuring 
line the absorption profile we can derive atmospheric tem- 
peratures, pressures and atomic sodium abundances. Those 
doublet profiles that have been measured show very differ- 
ent characteristics. The sodium absorption line profile in 
HD 209458b shows a narrow absorption core and broader 
plateau-like shape rather than smoot hly decaying wings 



Sing et al.ll2008al lbh. As explained by IVidal-Madiar et al.l 



2011bl ). the absorption depth should be greater towards 



smaller bandwidths due to an increase of the absorption 
cross-section. A plateau region could indicate an abundance 
drop where the increase in the cross-section is compensated 
for by the decrease in atomic sodium abundance in the al- 
titude regions sensed by narrower bands. The observation 
of pressure-broadened line wings shows that the Na 1 line is 
observed deep in the atmosphere and suggests a relatively 
cloud or haze-free atmosphere. 

Further analysis of tr ansit o bservat ions of HD 209458b 
by IVidal-Madiar et ail (|2011bl ) and IVidal-Madiar et all 
|2011a ). combining archive Hubble Space Telescope (H ST) 
and ground-based measurements l|Snellen et al.l l2008h , al- 
lowed a temperature-pressure (T-P) profile to be con- 
structed from 10 -3 to 10~ 7 bar. The shape of the line ab- 
sorption is characteristic of the atmospheric temperature 
profile, as higher temperatures increase the scale he ig ht thu s 
causing a steeper line slope. IVidal-Madiar et al.l (|2011aT l 
found that the temperature rose with altitude over the range 
10" 3 -10" 7 bar to 3600 ± 1400 K at the highest altitudes, 
a sign of the planet's thermosphere. Models suggest that 
very hot thermospheres should be common in close-in ex- 
trasolar planets, due to the lar ge amounts of UV radiation 
receiv ed fro m the nearb y star. iLecavelier des Etangs et al.l 
(|2004T ) and lYellel { 20041 ) . modelling the upper atmosphere 
and atmospheric escape of HD 209458b, showed that tem- 
peratures in the the rmosphere could r each u p to 10,000 K. 
Later models from iTian et"aH l|2005h and iGarcfa Mufiozl 
(|2007h showed temperatures of a similar order of magnitude. 

In contrast to the measured profile for HD 209458b, 
WASP-17b appears to show o nly a narrow abso rption core, 
with no observable line wings l|Wood et al.ll201ll ). This indi- 
cates that there could be absorbing species very high in the 
atmosphere, obscuring the line wings. 

For HD 189733b, broad-band spectroscopy has shown 
a spectrum do minated by Rayle i gh scattering from 3000 
to 10500 A dPont et all 120081 ; ISing et all I2009I . 1201 ll ; 
iGibson et al.ll2012f ). with longer wavelengths probing deeper 
in the atmosphere. This featureless broad-band spectrum 
likely indicates that the Na I line wings produced by pres- 
sure broadening deep in the atmosphere are obscured by 
continuum absorption higher in the atmosphere. Na I has 
been detected in HD 189733b from the grou nd, with a rela- 
tive a bsorption depth of (67.2 ±20.7) x 10" 5 |Redfield et all 
2008) compared to the continuum. This detection was later 
revised to (52.6 ± 16.9) x 10" 5 by I Jensen et all (|201lT ). 

An upper atmospheric temperature-pressure profile has 
only been derived for one exoplanet, HD 209458b. The aim 
here is to conduct a similar study for HD 189733b. These 
two planets are in many respects at opposite ends of the 
spectrum of hot Jupiters. HD 209458b has a bloated radius, 



displays a stratospheric temperature inversion and orbits a 
star that is much less active than HD 189733. In contrast, 
HD 189733b does not have a bloated radius, shows no strato- 
spheric temperature inversion and orbits one of the most ac- 
tive extrasolar planet host stars. HD 189733 is a somewhat 
cooler star than HD 209458 (4980 vs 6075 K). Additionally, 
the transit absorption spectrum of HD 209458b is domi- 
nated by atomic absorption from alkali lines, whereas that 
of HD 1 89733b appears do minated by a Rayleigh scatter- 
ing haze. iHeng et al.l (|201ll ) demonstrate that the presence 
of clouds and hazes can have very different effects on the 
temperature-pressure profile depending on their thicknesses 
and opacities. IR absorption due to a haze layer can heat the 
lower atmosphere and cool the upper atmosphere, whereas 
optical scattering will have the opposite effect, warming the 
upper atmosphere. The many differences between the two 
planets makes HD 189733b an ideal planet to study in de- 
tail for comparison with HD 209458b. 

Since many factors can affect the energy budget of a 
planetary atmosphere, such as advection efficiencies and 
opacities in the atmosphere, it is difficult to say in ad- 
vance how incoming stellar radiation on HD 189733b af- 
fects the atmospheric temperature profile. Differing lev- 
els of radiation from the host stars may have a signifi- 
cant effect on the upper atmosphe ric chemistry as well. 
iKnutson. Howard &: Isaacson! l|2010T l suggest that planets 
that orbit the most active stars have non-inverted strato- 
spheres and planets that orbit quiet stars do have inverted 
stratospheres. The authors suggest that the excess EUV ra- 
diation from a more active star can break down the com- 
pounds responsible for stratospheric temperature inversions, 
thus altering the T-P profile significantly even at altitudes 
lower than the regime observed here. 

In this paper, we present HST Space Telescope Imaging 
Spectrograph (STIS) medium-resolution transmission spec- 
tra of HD 189733b with the G750M grating. The increased 
resolution and high signal-to-noise of these spectra com- 
pared to the ground-based detection allows us to measure 
the spectral absorption depth profile of the sodium doublet 
and construct a vertical temperature profile of the atmo- 
sphere. Making an assumption of the pressure at a reference 
altitude, this translates into a temperature-pressure profile 
which can be compared to theoretical predictions. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Light Curve Fitting 

We observed 7 primary transits of HD 189733b using the 
HST STIS G750M grating (HST GO 11572, PI: David 
Sing). Three of these transits obtained data of high qual- 
ity. Of the others, two suffered from a positioning error of 
the spect rum when us i ng su b-arrays on the CCD, as also 
noted in I Brown et al.l |200ll '). Two further visits obtained 
no data as those visits suffered from guide star problems. 
Table Q] gives details of the 7 visits. 

Each dataset consists of a series of 144 spectra, each 
covering a wavelength range of 5808-6380 A. The G750M 
grating has a resolving power of R — 5000, which gives a 
resolution of ~ 1.2 A or ~ 2 pixels at 5890 A. The obser- 
vations were made with a slit of 2" width to minimise light 
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Visit 


JD 


Notes 


1 


2455135.31 


Positioning error. Spectrum blueward of 






6000 A used in part of the analysis 


2 


2455141.42 


No data acquired 


3 


2455148.62 


Whole spectrum useable 


1 


2455184.11 


No data acquired 


5 


2455423.72 


Positioning error. Spectra were 






unuseable 


6 


2455470.31 


Whole spectrum useable 


7 


2455530.22 


Whole spectrum useable 



Table 1. Table of HST STIS visits, showing data quality. 



Visit 3 



Visit 6 



Visit 7 



Fit Position 


no 


yes 


no 


yes 


no 


yes 


Parameters? 














DOF 


69 


66 


79 


76 


86 


83 


NFree 


8 


11 


8 


11 


8 


11 


"K.v 


1.948 


1.751 


5.705 


1.958 


3.632 


1.936 


BIC 


349 


167 


447 


183 


265 


211 



Table 2. Table of fitting statistics depending on whether the 
remaining x and y offsets were included in the fit. 'NFree' is the 
number of free parameters. It can be seen that there are significant 
improvements for all the visits when these parameters are fitted. 



loss. The data were then bias-subtracted, dark-subtracted 
and flat-fielded using the CALSTIS pipeline. The spectra 
were extracted afterwards using IRAF with an aperture 13 
pixels wide. IRAF was also used to remove cosmic ray con- 
tamination and to measure the spectral trace (the position of 
the spectrum on the CCD). These spectra were then shifted 
to a common wavelength scale by cross-correlating with the 
mean of the shifted spectra, and were then corrected to a 
rest frame by comparison with a model spectrum. 

Transit light curves were produced by summing the 
photon flux over the full spectrum in each exposure (see 
Figure [Q , and were mode l led u sing the analytical transit 
models of iMandel fc Agoll i|2002T) . Th e light curves exhibi t 
all the instrumental effects noted by [Brown et al.l i|200ll ). 
The systematic effects were corrected for by discarding the 
first orbit of each visit and the first exposure of the retained 
orbits, and then removing a fourth-order polynomial depen- 
dence of the fluxes with HST orbital phase based on a fit 
to the out-of-transit exposures. The Levenberg-Markwardt 
least-squares tec hnique was used for the fit using the IDL 
MPFIT package (|Markwardtj 2009') using the unbinned data. 
The in-transit sections obviously crossing stellar spots were 
removed in these fits. 

We also investigated whether fitting for any remaining 
x and y offsets with of the spectra on the detector (such as 
due to target motion within the aperture) with linear func- 
tions were justifi ed by improvement in the fit as done by 
ISing et all (|201ll ). These variations improved the reduced 
X 2 and significantly improved the Bayesian information cri- 
terion (BIC), which penalises models with larger numbers 
of free parameters, and prevents over-fitting the data. Table 
Ogives details of these improvements. 

Fitting the systematic trends in this way produced S/N 
values for the white light curves between 10,000-11,000 (pre- 
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Figure 1. White light curves for (top to bottom) visit 3, visit 
6 and visit 7 with the systematic effects removed. An arbitrary 
flux offset has been applie d for clarity. For each v isit, the best fit 
analytical transit model of Mandel & Asol (2002) is over-plotted. 
Plotted underneath are the residuals for (top to bottom) visit 
3, visit 6 and visit 7 after the systematics have been removed. 
An arbitrary flux offset has been applied for clarity. The regions 
where the planet crosses stellar spots can be clearly seen. 



cision levels of 90-100 ppm), which are 78-86 per cent of 
the Poisson-limited value. These per exposure precision lev- 
els are an improvement on the results when not fitting for 
these additional position-dependent parameters (with preci- 
sion levels of 130-160 ppm). 

We checked for red noise by observing whether residuals 
binned in time followed a N' 1 ^ 2 relationship, where TV is the 
number of points in the time bin. We found no significant 
evidence of red noise in visit 3 or visit 6 for the white light 
curve. For visit 7, the magnitude was 3 x 10~ 5 for the white 
light curve, which corresponds to a S/N of ~ 33,000. It 
translates into an error in depth measurement in the planet's 
atmosphere of ~ 50 km, which is less than 1/4 of a scale 
height in the lowest regions (smallest scale height) probed 
by these observations. Since our uncertainties are ~ 100 — 
300 km in altitude, these low levels of systematic trends 
do not dominate the uncertainties in our measured spectral 
absorption depth profile. Smaller bandwidths are generally 
closer to Poisson-limited, and it was also found that the red 
noise m a 20 A band was negligible (< 1 x 10 J ). 

We also used linear regression to check for remaining 
correlation between fitted parameters, and all correlations 
were of very small order (in the range 0.0008-0.0001). Fig- 
ure[T]also shows the residuals after removing the polynomial 
characteristic of the systematic effects. Visits 3 and 6 still 
exhibit structure in the transit due to occulted starspots, 
where the flux is higher than the transit model predicts (see 
Section 13. 3[) . It was observed that the transit depth and 
hence the sodium absorption depths for visit 6 was shal- 
lower than the other two visits. This could be due to an 
occulted starspot in the centre of the transit in addition to 
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those at the edges, so that the planet would cover part or 
all of at least one starspot during the whole transit. 

2.2 Stellar Limb Darkening 

Correcting for the effects of stellar limb darkening on the 
transit light curves and the spectral absorption profiles is 
important at visible wavelengths. We compared two limb 
darkening corrections based on two stelar models; the Ku- 
rucz (1993) ID ATLAS stellar atmospheric model 

and a 

fully 3 D time-dependent hydrodynamic stell ar atmospheric 
model jSing et al.l l201ll ; lHavek et all I2012T ). The models 
were compared using visit 7, where there seem to be no 

occulted spots at t he limbs of the tr ansit. 

In the ID case. lSing et~aH (|2009l ) and lSind (|201Ch found 
that a four-parameter law broke down for small [i, where 
H — cos(#), and 9 is the angle in radial direction from the 
disc centre. Since the first coefficient has a greater effect 
at small angles, it was not used, and we instead fitted a 3 
parameter law: 

*M = 1 - c 2 (l - u) - ca(l - m 3/2 ) - C4(l - m 2 )- (1) 

The stellar parameters used were T e g = 5000 K, log(g) = 
4.5, and [Fe/H]= 0.0, along with the transmission function 
of the G750M grating. The resulting coefficients for white 
light were c 2 = 0.8918, c 3 = 0.0515 and c 4 = -0.1930. This 
fit gave a reduced \ 2 of 2.52 and BIC of 268. 

In the 3D case, the stellar limb darkening was fitted 
with a 4-parameter law, as it appears to perform well at 
small fj,: 

= l- Cl (l- M 1 / 2 )- C2 (l- M )- C3 (l- Al 3 / 2 )-c 4 (l-/i 2 ).(2) 

The parameters for white light were ci = 0.7043, c 2 = 
-0.4493, c 3 = 1.0538, and c 4 = -0.4569. It was found that 
the fits and the magnitudes of red noise improved over the 
ID case. The reduced \ 2 when using the 3D model was 1.94, 
with a BIC of 211. Despite this, however, the measured 
planet -to-star radius changed by less than la. ISing et alJ 
(j201lT ) give an in-depth discussion of the comparison be- 
tween the ID and 3D cases for broadband G430L transit 
observations (2900-5700 A). They also find that the fits to 
their data are significantly better using the 3D model than 
the ID model. 



3 ANALYSIS 

3.1 Spectral Absorption Depth Profile of the Na I 
Doublet 

We measured the spectral absorption depth profile ver- 
sus wavelength around the sodium lines by comparing 
the absorption depths in the sodium wavelengths to ab- 
sorption depths in a reference band during transit. For 
the reference band , we u sed the 'wide' band defined by 
ICharbonneau et al.l l|2002r i , which is a combination of a blue 
(5818-5843 A) and a red (5943-5968 A) band. The three 
complete visits were used for this analysis, as was visit 1, 

1 See http://kurucz.harvard.edu/stars/hdl89733 



since we only required wavelengths up to 5968 A. We did not 
use the fitted radius from visit 1, only the differential flux 
measurement described here, so this visit is not discussed in 
Section 12.11 We tested using the 'medium' reference band 
as defined by Charbonneau et al. (2002) as well as a refer- 
ence band closer to the doublet, and found that the choice 
of reference band shifts the profile in altitude by less than 

This differential method has the advantage that system- 
atic errors largely cancel out, as they are largely wavelength 
independent, giving a robust way to measure Na I absorp- 
tion. It also averages the blue and red background compo- 
nents, giving an average that largely compensates for the 
effects of limb darkening in the sodium bands. This method 
is also particularly useful here due to the presence of oc- 
culted starspots at the limbs of two of the visits, as the 
flux variations due to transiting in front of occulted spots 
are similar in the Na band and the 'wide' bands, and thus 
cancel out when comparing the bands. 

However, this method does not account for differential 
limb darkening in the sodium line core wavelengths. The 
limb darkening is weaker at the line core compared to the 
continuum, because the core of the stellar Na I line is pro- 
duced above the photosphere, further away from the cen- 
tre of the star than where the continuum is produced. At 
these altitudes, the temperature gradient becomes smaller 
(see Hayek et al., submitted). The differential limb darken- 
ing leads to a slight underestimation of the absorption depth 
in the line core, since the limb darkening effects are not fully 
cancelled out by the differential method. This effect was cor- 
rected by comparing limb darkening models for each of the 
evaluated sodium wavelengths with that of the 'wide' band. 
These models were given a fixed planetary radius, to eval- 
uate the difference in model limb darkening, and then the 
differences were added to the sodium absorption depths. The 
effects were found to be negligible for wavelengths greater 
than 3 A away from the line cores. 

Figure f2] shows the spectral absorption depth profile 
around the Na I feature binned to the instrument reso- 
lution (2 pixels) and Table [3] lists the differential absorp- 
tion depths for each wavelength bin. We find in Section f3. 31 
that the effects of starspots on this profile are minimal and 
within our observational uncertainties. Comparin g Figure [2] 
to the e quivalent spectrum for HD 209458b from Sin g et alJ 
(|2008af ) shows that the measured differential absorption for 
HD 189733b is deeper in the line cores. There is one data 
point in the centre of the two doublet lines which appears to 
have an unexpectedly low absorption in the unbinned data, 
and can still be seen in the binned data. It is not associ- 
ated with any excess red noise, and so we could not scale 
the uncertainty any higher. This a nomaly was also ob served 
in the HD 209458b spectrum from lSing etafl (|2008Vj ) using 
the same instrument, but there are no obvious defects on 
the CCD that could account for this effect. 

3.2 Integrated Absorption Depth Profile of the 
Na I Doublet 

In order to allow direct comparison with the previously 
published Na I p r ofiles of HP 2094 58b and WASP-17b 
l|Sing et all l2008al : IWood et ail l201ll ). we also compute 
the differential integrated absorption depth profile as a 
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19.04 


5871.61 
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21.41 
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21.21 
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-34.05 
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5876 04 
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5877.15 


-17.25 


20.63 
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6.80 


21.09 
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31.16 


23.01 


5880 47 


34.76 


20.70 
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-8.91 


21.23 
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24.54 
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21.49 
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30.13 
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20.34 
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18.51 


21.61 
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57.24 


24.06 


5889 33 


40.69 


36.82 




213.89 


36.17 


5891 55 


86.23 


24.62 


5892 fifi 


31.31 


26.42 


5893 77 


-7.95 


23.81 


58Q4 87 


48.48 


29.33 


5895 98 


141.22 


40.31 


5897.09 


59.56 


25.44 


5898.20 


71.53 


20.53 


5899.31 


-21.23 


25.01 


5900.41 


5.60 


21.23 


5901.52 


41.71 


21.60 


5902.63 


-2.19 


23.29 


5903.74 


10.17 


21.49 


5904.85 


25.18 


22.27 


5905.95 


7.36 


21.65 


5907.06 


43.25 


20.97 


5908.17 


-10.27 


21.61 


5909.28 


-65.38 


18.88 


5910.39 


6.79 


23.31 


5911.49 


-27.29 


19.54 


5912.60 


10.12 


19.95 


5913.71 


13.50 


22.14 


5914.82 


-2.03 


26.01 


5915.93 


-11.71 


23.66 


5917.03 


-46.86 


21.42 


5918.14 


12.27 


19.40 


5919.25 


24.25 


22.80 









Table 3. Weighted average of spectral absorption depth (AD) 
in each wavelength (A) around the Na I doublet compared to 
absorption in adjacent bands. The data are binned to the STIS 
instrument resolution of 2 pixels. 
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Figure 2. Spectral absorption depth profile at each wavelength 
in the sodium region with la error bars, binned to the STIS reso- 
lution of 2 pixels. The Na I Dl and D2 doublet lines are resolved, 
as shown by the largest two peaks. The absorption depths for the 
doublet are greate r than the equivalent profile for HD 209458b 
ijSing et al.ll2008ah . 



function of bandw i dth. U sing; the method described by 
ICharbonneau et al.l (|2002t ). the absorption depth during 
transit in a band centred on the Na I doublet was com- 
pared the same blue and red reference bands as used for 
the spectral absorption depth profile. The bands centred on 
each sodium doublet line (5896 and 5890 A) were then in- 
creased in size from 3 to 80 A, with a single band being used 
where the bands were large enough to encompass both lines 
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Figure 3. The integrated absorption depth profile comparing the 
absorption depth in the Na I doublet to the absorption depth in 
adjacent bands as a functi on of bandwidth ar ound the sodium 
doublet. The detection from I Jensen et al] <|201 lT > is shown in blue 
with a filled square. Our more precise results agree with their 
measurement at the 1 a level. 



together (bandwidth ^ 12 A). We heareafter refer to this 
profile as the integrated absorption depth profile. 

Figure [3] shows the integrated absorption depth pro- 
file as a weighted average for the 4 visits. The line appears 
wider than in reality due to a dilution effect of increasing 
the bandwidth, weakening the sodium absorption signature 
confined to smaller bandwidths. The profile gives a depth of 
(51.1 ± 5.9) x 10" 5 for a bandwidth of 12 A, which was the 
bandw idth used by iRedfield et all (|2008l ) and J ensen et al,l 
(2011) for their ground-based measurements. Our result is 
within 1 a of these previous results, confirming the absorp- 
tion depth found with an increased S/N of almost a factor 
of 3. 



3.3 The Effect of Occulted Starspots 

The local flux increases seen as bumps in the light curves 
in Figure [1] are most likely due to occulted starspots. 
HD 189733 is a very active star, showing periodic flux varia- 
tions over time, and star spots are clearly s een in other tran- 
sit light curves as w ell (|Sing et all l201ll ; iPont et al.l 120081 ; 
iHenrv fc Winnll2008h . Since starspots are darker than the 
surrounding stellar surface, the planet blocks less of the stel- 
lar light when crossing a spot than when crossing a non- 
spotted region of the surface, resulting in a local flux in- 
crease in the light curve compared to the expected model. 
Occulted spots could therefore have an effect on the mea- 
sured sodium absorption depth because they decrease the 
inferred planetary radius. 

The differential method of measuring the absorption 
profile should mostly cancel out the effects of stellar spots. 
However, we do not a priori know how much the level of 
unocculted stellar spots on the surface affects the stellar 
sodium signature, and hence the effect on the inferred plan- 
etary absolute absorption depth. This effect could be greater 
in the line core than in the broad regions of the line and could 
be significant if the sodium signature is strongly dependent 
on the spot temperature. The sodium data here are used 
to empirically determine the effect of occulted spots on the 
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Figure 4. Weighted average of integrated absorption depth in 
a band centred on the Na I doublet compared to absorption in 
adjacent bands. Black: using all exposures, blue: using only expo- 
sures clearly not contaminated by occulted starspots, red: using 
only exposures that appear to contain occulted starspots. This 
shows that the starspots do not have a significant effect on the 
Na 1 absorption profile measurement at our level of precision. A 
colour version is available in the online journal. 



differential Na I spectral absorption depth profile. We per- 
form this analysis on the integrated absorption depth profile 
because the effects will be more severe in this profile than in 
the spectral absorption depth profile since they will be most 
pronounced in the line core. Figure [4] shows the mean inte- 
grated absorption depth versus bandwidth for the 4 visits 
compared with the mean integrated absorption depths eval- 
uated using only the exposures that appear to occult stellar 
spots, and also only those that do not. 

As expected, the derived integrated absorption depths 
are smaller for when evaluated from exposures occulting 
spots, with a value of (39.7 ± 8.2) x 10 -5 integrated over 
a 12 A band, compared to (56.3 ± 9.6) x 10 -5 derived from 
the exposures not occulting spots. The effect at even smaller 
bandwidths, where we expect the difference to be more no- 
ticeable, is around 5— 15 x 10 -5 , still within our observational 
error bars. For our purposes, we can therefore assume that 
the spots have a similar spectrum to the non-spotted sur- 
face. In this case, unocculted starspots are expected to have 
an even less significant effect than occulted starspots, so we 
can assume that these effects are also smaller than our uncer- 
tainties in absorption depth. It has, however, been noticed 
that at both longer and shorter wavelengths in the broad- 
band data, the effects become more significant (S ing et "all 

in!]). 



3.4 Broadband 5808-6300 A Spectrum 

The G750M data can also be used to look for other spectral 
features and to determine the shape of the overall spectrum 
from 5808 to 6380 A. Spectra were created by binning the 
spectral time series into 5 wavelength bins of ~ 110 A for 
visits 3, 6 and 7. Each band was fitted individually in the 
manner described for the white light curve. The fitted planet 
to star radius ratios were then used to construct a binned 
transmission spectrum. 

Determining accurate radii required the light curves to 
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Figure 5. Transit exposures sampling stellar spots ovcrplottcd 
with the best fit spot model for visits 3 (top) and 6 (bottom). 
A non-spotted model has been subtracted from both. A colour 
version is available in the online journal. 



be corrected for occulted stellar spots. The fluxes of each 
exposure contaminated by occulted spots in visits 3 and 6 
were measured with the white light curve, producing the 
shape of the features. These were then used to model the 
spots in each of the wavelength bands, by fitting the fea- 
tures in each band with the shape measured from the white 
light curve and an amplitude parameter. This assumes that 
the shape of the spots is roughly the same in all the differ- 
ent bands and that the spectra of the spots are constant. 
Figure [5] shows the spotted exposures overplotted with the 
best fit spot amplitudes for each wavelength band. Figure [6] 
shows the data residuals once the spots have been removed. 
Comparison with Figure [1] shows that the spots have been 
effectively removed. 

Starspots are cooler than the surrounding surface, and 
so are comparatively brighter at longer wavelengths. Oc- 
culted spot features are t herefore more sev ere at blue wave- 
lengths than at red fe.g. IPont et al.ll2007h . Figure [7] shows 
the best-fit spot amplitude parameters along with a pre- 
dicted slope for spots of temperature 4000 K. Sin g et all 
(|2011f ) found that this temperature best matched their data 
over a broad wavelength range. The amplitude was defined 
to be 1 for white light. Our slope is in agreement with this 
model, but cannot be used to constrain the spot tempera- 
ture further than the constraints provided by previous work. 
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Figure 6. Residuals of the data minus spotted models for visit 
3 (top) and visit 6 (bottom) for the G750M band divided into 5 
~ 110 A bins. Arbitrary flux offsets have been applied for clarity. 
The colours and wavelength bands are the same as in Figure [5] 
and a colour image can be found in the online journal. The flat- 
ness of the residuals show that the occulted starspots have been 
effectively removed. 
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Figure 7. Average best fit spot magnitudes from the models for 
visit 3 and visit 6, plotted against band. It can be seen that there 
is a wavelength-dependent slope, but this is within the errors of 
the measured parameters. Overplotted in red is a model spectral 
signature of occulted spot amplitude, assuming that the spots 
have a temperature of 4000 K. The measured slope is in agreement 
with this model, but cannot be used to constrain the temperature 
any further. 

It can be seen that the amplitude increases slightly towards 
blue wavelengths as expected, although this is only a 1 a 
effect over our wavelength range. 

The planetary radius and spot magnitude parameter 
were initially fitted together and then fitted one at a time 
iteratively until the solutions converged, allowing the radii 
in each wavelength band to be fitted while taking into ac- 
count the shape of the spots in the light curve. The spec- 
trum was then corrected for u nocculted spots b y reducing 
the measured radii by 1% (see ISing et alTl201ll ) and aver- 
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Figure 8. The STIS G750M broad spectrum overp lotted with 
a Rayleigh scattering prediction from the ACS data l|Pont et alj 
2008; Sing et ah 2011 ) at the temperat ure of 1340 K determined 
by iLecavelier des Etangs et al.l ll2008al) (solid line) with a range 
of ±150 K (dashed lines). 



Wavelength (A) 


Rp/R* 


Error 


5809-5922 


0.15638 


0.00025 


5923-6037 


0.15632 


0.00019 


6040-6150 


0.15618 


0.00034 


6150-6264 


0.15601 


0.00025 


6264-6378 


0.15612 


0.00016 



Table 4. Fitted planet-to-star radius ratios for ~ 110 A bands, 
used to construct the broad-band spectrum for the G750M band. 
The radii have been corrected for occulted and unocculted stellar 
spots. 

aged for the 3 visits. Figure [8] shows the resulting spec- 
trum and Table [4] shows the fitted values of Rp / R* . Fig- 
ure [5] also shows a Rayleigh slope for an atmospheric tem- 
perature of_1340^J^0_J£_dete£mined^^m previous obser- 
vations (jLecavelier des Etangs et al.l 12008a). Our data are 
consistent with Rayleigh type scattering, but we cannot fur- 
ther constrain the uncertainty on the measured temperature 
with our data, although we do reach precisions of approxi- 
mately one atmospheric scale height for each of our binned 
spectral points. Additionally, we see a flat spectrum with 
spectral bins narrower than those of either of the two previ- 
ously observed broad-band spectra, which is consistent with 
the conclusion that the Na I line wings are absent. It also 
indicates that there are no other broad spectral features in 
our wavelength range. 

In contrast to the flat spectrum measured here, the 
spectrum of HD 209458b showed significant excess a bsorp- 
tion at 6200 A (|Sing et al.ll2008al ; I Desert et ai]|2008D . The 
non-detection of significant excess absorption at 6200 A in 
HD 189733b suggests that, either this feature is real or any 
systematics responsible for the feature vary with time. 

3.5 Determining Upper Atmospheric 
Temperatures 

Visible transmission spectroscopy and the sodium D doublet 
in particular should probe the upper atmosphere, at pres- 
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sures below 150 mbar to as low as 10 9 bar, where we could 
see very high t e mperatures ( Lecavelier des Etangs et alj 
20041: lYelld 12004 iGarcfa Munozl 120071 : iMoses et al.ll201ll ). 
Vidal-Madiar et alj (|201 lbl laf) measured a temperature of 
3600 ± 1400 K for the upper atmosphere of HD 209458b, 
using the narrowest bands (highest altitudes) of the exo- 
planet's sodium absorption profile. Here, we use the same 
method to determine the temperature ranges covered by our 
measurements. We generate models of the sodium D doublet 
absorption lines at various temperatures and calculate spec- 
tral absorption depth profiles at the spectral wavelengths 
of the measurements. This method is independent of abso- 
lute absorption depth and depends only on the absorption 
profile slope, which means that we do not have to know a 
reference pressure scale or the abundances to measure lo- 
cal temperatures. We assume constant sodium abundance, 
which is valid for local measurements. 

The atmospheric altitude, z(A), that can be ascribed 
to the observed absorption as a function of wavelength in 
a transmission spectrum, is given by the equation from 
ILecavelier des Etangs et ah! (|2008aT ): 



a(A) = 



z(X) 



kT, 



fig I r cq ^ kTfig J 



(3) 



where k is Boltzmann's constant, T is the temperature, /x 
is the mean molecular weight of the atmospheric compo- 
sition, g is the surface gravity, £ is the elemental abun- 
dance, <t(A) is the wavelength dependent absorption cross- 
section, P is the pressure at the reference zero altitude, 
r cq is the optical depth at the transit radius, shown to 
be approximately constant by Lecavelier des Etangs et al. 
(2008), and _Rp is the radius of the planet. This relation- 
ship assumes hydrostatic equilibriu m and the ideal gas law. 
ILecavelier des Etangs et alj (|2008aT ) show that the temper- 
ature can be inferred from the slope of the line by differen- 
tiating Equation @: 



dz_ _ kT <9(ma(A)) 
dX ~ 



I if) 



dX 



M g f dQna(X)) 
k V d\ 



— ■ (4) 



In the case of the sodium doublet absorption, measuring the 
slope of the spectral absorption depth profile in the actual 
data can be used to derive temperatures for different parts 
of the profile, and thus temperatures for different altitudes 
where the absorption takes place: deeper in the atmosphere 
from the broader parts of the lines an d at higher altitudes 
from t he narrower parts of the lines fe.g. lVidal-Madiar et al.l 
l2011bh . Any non-constant slope of the measured spectral 
absorption depth profile indicates that the profile probes 
more than one temperature regime. 

For HD 189733b, we assume g — 2141 cms - and an at- 
mosphere composed of mainly hydrogen and helium (85 per 
cent H, 15 per cent He), where /i = 2.3 x 1.6726 x 10~ 24 g. 
The sodium abundance, £, planetary radius, Rp, and the 
reference pressure, P , do not affect the slope of the model 
spectral absorption profiles and hence the temperatures. 
We arbitrarily set t hese values to ^Na/^H = 1.995 x 10~ 6 , 
a solar abundan ce jLoddersl 2003 | ). Rp = 1 . 138 R j, and 
P = 150 mbar. ILecavelier des Etangs et al.l l|2008al ) show 
that T cq = 0.56. To determine <r(A) for the sodium doublet, 
each line was modelled as a Voigt profile, where 



H = 



/ 



H, 



(5) 



and Oo is the absorption-cross section at the centre of the 
relevant line, Aud is the Doppler width, H is the Voigt 
function, which includes functions for natural, Doppler and 
pressure broadening, e is the electronic charge, m e is the 
mass of the electron, and / is the absorption oscillator 
strength. We used / = 0.6405 for the D 2 line and / = 0.3199 
for the Di line (Steck, 20100. The wavelength-dependent 
cross-section, <r(A), was calculated separately for the Di 
and D2 lines, and the combined profile was worked out us- 
ing <t(A) = (t(A)di + o"(A)d2- The Doppler width is given 
by Auo = Vo/cy // 2kT/ fiNa, where v„ is the central fre- 
quency of the doublet being computed, c is the speed of 
light and fiNa is the mean molecular weight of sodium, 
23 x 1.6726 x 10 -24 g. We used the central wavelengths from 
Steck (2010), with X OD1 = 5895.98 A and X OD2 = 5890.00 A. 
These also match the central wavelengths observed in the 
data. 

We used Equation (|3]) to calculate a spectral absorption 
depth profiles using er(A) and fitted these to each region of 
the lines with T as a free parameter. To allow us to fit each 
region locally and account for global changes in abundance, 
we also fitted the profiles with a free altitude parameter, z , 
which allows the profile to shift up and down. 

We find two temperature regions. Fitting the whole 
Na profile with one temperature gives a \ 2 °f 55.5 and 
T = 2000 K for 35 DOF with a BIC of 59.0. However, fitting 
for two temperature regimes, one at low altitudes probed by 
the 'line wings' and one at high altitudes probed by the 'line 
cores' results in a superior fit, with a BIC of 44.6. The sig- 
nal in the 'line wings' (5874-5886 A and 5899-5012 A) of the 
Na I D lines is significant at the 2.9 sigma level added over 
the whole wavelength region for all visits, indicating a likely 
detection separate from the 'line cores ' (5887-5898 A). The 
best fits are shown in Figure [9] with the heights in km be- 
ing relative to the continuum and showing the atmospheric 
regions probed by each part of the spectral lines. 

The best fit to the 'line wings' of the Na I D feature is 
1280±240 K, with \ 2 = 28.0 (24 DOF). This region approx- 
imately covers altitudes of less than 500 km above the con- 
tinuum, which is within the same vertical region in the at- 
mosphere as the b lueward rise in ab sorption depth observed 
in the near-UV bv lSing et al.l (|201ll ). who found a tempera- 
ture of 2100 ±500 K. This marginally suggests that another 
mechanism may be required to explain the excess absorption 
in the blue region of the spectrum, though the temperatures 
agree at the 1.5a level. The best fitting temperature for the 
Na I 'line cores' was found to be 2800±400 K, with x 2 = 11.0 
for 11 DOF. This region corresponds approximately to alti- 
tudes of 500-4000 km above the continuum. 

We fit the data neglecting pressure broadening. Includ- 
ing pressure broadening in the fit with pressures of 60- 
150 mbar for the line wings (see Section 14.2. 1[) has very 
little effect, with the best fitting temperature increasing by 
~ 20 K and the \ 2 increasing to 32.6 for 24 DOF. Pressure 
broadening does not affect the innermost core of the line. 

In order to consider the spectral dispersion in the fits, 
we moved the centres of our fitted sodium model by one 

2 See http://steck.us/alkalidata 
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Figure 10. Temperature-altitude profile for HD 189733b. The 
altitude values shown are binned to the HST STIS resolution, 
relative to the continuum level. The thick vertical lines indicate 
isothermal atmospheric layers. The horizontal error bars show the 
range of temperatures that match the spectral absorption depth 
profile for the given atmospheric region. The lowermost point is 
the temperature derived from the Rayleigh continuum slope by 
ILecavel ier des Etangs ~et alj l2008al ). The temperature rise with 
increasing altitude that we observe occurs much higher than the 
stratosphere. 
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Figure 11. Spectrum over the G750M band (black squares) com- 
pared with previous work. Also shown are points from the G750M 
spectrum at 2, 5, 12 and 18 A from the centre of the doublet 
(black circles) . The red l ines show the Rayleigh scattering p rcdic- 
tion from the ACS data (Lccavclicr des Etangs ct al. 2008a). The 
blue line shows a m odel assuming a haze-free atmosphere from 
iFortnev et al] §010) at solar Na I a bundance, normal ised to the 
radius at infrared wavelengths from Agol ct al. J2010f) ■ Blue cir- 
cles show model values binned to the G430L and ACS resolution. 
A colour version is available in the online journal. 



pixel (0.54 A) bluewards and redwards of the centre val- 
ues and found that the 'line core' temperatures changed by 
only ~ 200 K. We also calculated the Doppler broadening of 
the line that should result from the planet's motion about 
the star, which could lower the temperatures derived from 
the line shape. We found that, for the phase of the observa- 
tions, the maximum blueshift (during ingress) and redshift 
(during egress) is < 0.7 A, which is below our instrumental 
resolution. 

We further tested the effect of adding Rayleigh scatter- 
ing at the lowest altitudes by combining the cross-sections 
for the sodium doublet and a Rayleigh signal. Keeping the 
'line core' temperature fixed at 2800 K and fitting the 'line 
wings' with a sodium and Rayleigh signal gives a \ 2 °f 
35.6 instead of 28.0 in the sodium only case, indicating 
that a Rayleigh component to the signal is likely negligi- 
ble . Since the temperature derived from the Rayleigh slope 
by iLecavelier des Etangs et alj |2008a) is the same as the 
fit obtained here by assuming that the whole signal in the 
wings is due to sodium, the presence of a Rayleigh compo- 
nent in our signal does not affect our calculated atmospheric 
T-P profile. 

The resulting temperature-altitude profile is shown in 
Figure 1101 where the altitudes correspond to the approxi- 
mate atmospheric height ranges probed by each of the fitted 
wavelength regions of the doublet. This profile includes the 
previously measured temperature o f 1340 ± 150 K found by 
ILecavelier des Etangs et al.l (|2008al ) for the Rayleigh con- 
tinuum slope in the sodium region, at altitudes of z = 0- 
200 km. 



4 DISCUSSION 

4.1 Interpretation of the Narrow Absorption 
Lines 

Figure 1111 shows th e prev i ously observed broad-band 
spectrum |Pont et al. 200$ ISing et al.l l201ll) compared 
to a model spectru m assuming a haze-free atmosphere 
|Fortnev et al.lr2010t ) and a solar Na I abundance. The ob- 
served spectrum lacks the broad alkali line wings and H2O 
signatures. The broadband data are binned to ~ 500 A, 
showing that the equivalent widths of those bins are not 
significantly higher than the featureless spectrum. Figure [TT] 
also shows our measured spectrum for the G750M band as 
well as points from the medium resolution spectrum at 2, 5, 
12, and 18 A from the centre of the doublet, confirming a 
narrow Na I feature. The previous broad-band spectroscopy 
would not have had sufficient spectral resolution to resolve 
such a narrow feature unambigu ously. The 3 ob servations 
from STIS G750M, STI S G430L (|Sing et al.ll2o7lh and ACS 
HRC (|Pont et al.ll2008T ) all independently agree that there 
are no observable broad Na I line wings. 

The species responsible for the Rayleigh scattering 
signature observed in HD 189733b is unknown, but the 
signature should place constraints on the composition of 
the scatte r ing s pecies (|Lecavelier des Etangs et ail [2008a; 
ISing et all l201ll ). Rayleigh scattering dominates only for 
particle sizes much smaller than the wavelength and where 
the scattering efficiency is dominant over the extinction ef- 
ficiency. H2, which is abundant in the atmosphere, is much 
smaller than the observed wavelengths. A narrow Na D dou- 
blet profile can be explained by a low sodium abundance 
which can hide the Na I line wings beneath the H2 Rayleigh 
signature. 

Alternatively, if the abundance is solar or above, a 
high-altitude material could cover the lower regions of 
the line, allowing only the upper regions to be observed, 
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Figure 9. Plots showing isothermal model fits to the different wavelength regions of the data. The best fitting models are shown in red, 
binned to the instrument resolution. Left: Fit to 'line wings ' (5874-5886 A and 5899-5912 A), which probes atmospheric regions less than 
~ 500 km above the reference altitude. The best fitting temperature is 1280 ± 240 K, a temperature that is similar to the temperature 
obtained by fitting the broad-band continuum absorption with Rayleigh scattering at the same wavelength range. Right: Fit to the 'line 
cores ', at 5887-5898 A. This wavelength region probes higher atmospheric regions, greater than ~ 500 km above the continuum. The 
temperature is found to increase to 2800 ± 500 K. 



where pressures are low and there is no observable pres- 
sure broadening of the feature. If the Rayleigh scatter- 
ing is from a haze of condensate particles, then the con- 
densate must be transparent enough for the absorption to 
be negligible relative to scattering, which rules out many 
condensates such as MgSiCU and o xygen-deficient silicates 
l|Lecavelier des Etangs et al.ll2008ah . One possible species is 
MgSiOs. Assuming a uniform haze layer with a single size 
component, the observed broad-band data from STIS and 
NICMOS constrain the particle size to be in the range 0.01- 
0.03 um by requiring that the molecule scatters at all ob- 
served wavelengths. 

The narrow range of particle sizes needed to be consis- 
tent with the observations shows that it is potentially diffi- 
cult to explain the broad-band Rayleigh signature over the 
4000-18700 A range with scattering by a haze of condensates 
of constant grain size. It is more likely that a condensate 
would vary in grain size with depth, and that there could be 
more than one absorber (e.g. aerosols) but since the scatter- 
ing species is so unknown, we present only a representative 
T-P profile based on these grain sizes for a single absorber. 



4.2 Temperature-Pressure Profile 

The temperatures of the upper atmosphere probed by the 
optical observations increase with increasing altitude. For 
our scale heights, we assume that most of the hydro gen is 
in the form of H 2 . At higher altitudes, lYelld (|2004h sug- 
gest that the scale height could increase further due to the 
change in dominance of hydrogen from Hb to H, which would 
decrease the mean molecular weight, fj,. A lower mean molec- 
ular weight would extend the line witho ut having to increase 
the temperature. iGarcia Munozl (|2007l ) show that H domi- 
nates over H2 only for pressures lower than P — 10 _s bar. 
iMoses et al.l (|201 lh also estimates that photodissociation of 
H2 occurs sharply at 10 -8 bar, although this pressure de- 
pends on where in the atmosphere the stellar radiation is 
deposited and could be at higher pressures. Our Na I line 



cores (highest altitudes) occur at ~ r/Rp ~ 1.06, where H2 
should dominate. 

Our measurements probe atmospheric regions well 
above the stratosphere, which could indicate that we are 
observing the base of the atmospheric thermosphere. To 
compare our temperature profile with theoretical models, 
we generate a T-P profile, using our derived temperatures 
to work out the corresponding atmospheric scale heights, 
and hence the pressure structure of the atmosphere relative 
to the reference level. We assume a well-mixed atmosphere 
where the scale heights derived from the sodium line will be 
representative of the atmosphere. 

The measured spectral absorption depths are only rel- 
ative to the absorption depth in the reference bands de- 
fined in Section 13.11 and hence the derived T-z profile 
can only be converted into a T-P profile if the pres- 
sure at the reference altitude is known (see discussion in 
iLecavelier des Etangs et' al. 2008b). Determining a reference 
pressure requires the species responsible for the continuum 
absorption in the reference bands and its abundance to be 
known. In the case of HD 209458b, the identification of 
Rayleigh sca ttering by H2 provided the required abundance 
information (|Vidal-Madiar et alj|20ilbh . Since the contin- 
uum absorption in HD 189733b is due to an unknown species 
of unknown abundance, the inferred pressure-altitude rela- 
tionship is not well constrained and can shift up or down 
depending on the reference pressure. 

We used a pressure profile based on a reference pres- 
sure of 150 mbar at z = (Rp = 0.15628) to determine 
a pressure-altitude relationship assuming H2 as the scat- 
tering species. Any other Rayleigh scattering species will 
be at higher altitudes, meaning that the pressure profile 
based on H2 as the scattering species is an upper limit. 
The T-P profile was arbitrarily shifted to lower pressures 
to consider T-P profi les that result from differen t possi - 
ble scattering species. ILecavelier des Etangs et all (|2008al) 
show that a haze of MgSiOs grains of sizes of 0.01-0.03 um 
at a temperature of 1340 K will be at pressures between 
10 -3 < P < 10~ 5 bar for solar magnesium composition. 



T-P Profile for HD 189733b 11 



4.2.1 The Terminator T-P Profile of HD 189733b 

Figure [121 shows the T-P profile derived assuming that the 
Rayleigh signature is due to a high-altitude silicate haze. 
The reference pressure is 10 -4 bar, representative of the 
10" 3 -10" 5 bar range that fits the broad-band observations. 
Depending on the species, this profile could shift signifi- 
cantly. In the case assuming instead that H2 is the scattering 
species, the derived pressures would be 150 mbar for the ref- 
erence level, 65 to 6.6 mbar for the lower region of the observ- 
able line at 1280 ±240 K and 6.5 to 0.02 mbar for the upper 
region of the line at 2800±400 K. We also show the tempera- 
tures and pr essures observed for lower altitudes (higher pres - 
sures) from IKnutson et al.l ll2007l); iKnutson et al.l (l2009h . 
ICharbonneau et al.l |200St ) and iDeming et al.l ^OOrj ). The 
temperatures in the last tw o works we re adju sted to the 
terminator temperature by H eng et al.l (|201 lh and these 
are the values we use. We also c ombi ne the tempera- 
ture measurement from Sing et al.l (120091 ^ with that from 
iLecavelier des Etangs et al.l ( 2008al ) since they are both de- 
rived from fits to the broad-band Rayleigh slo pe. The models 
shown for compariso n to HD 189733b are from Fortnev et al.l 
l|201Q] ). lYellel <|2004h and iGarci'a Mufiozl j2007h . Both upper 
atmospheric models were computed for HD 209458b. The 
vertical bars which indicate pressure range in the figure for 
HD 189733b do not include the uncertainty in the pressure 
reference level. They indicate the altitude range over which 
we match each temperature. 

The pressures where we observe the strong temperature 
rise are at least an order of magnitude higher than the pres- 
sures of the base of the thermosphere found in the models 
of the upper atmosphere of HD 209458b, assuming a high- 
altitude haze as a reference level. This difference increases 
even further, by 3 orders of magnitude, if we assume H2 
as the scattering species. The discrepancy between our data 
and the models could indicate that we are observing a lower- 
altitude temperature rise rather than the temperature rise 
characteristic of the base of the thermosphere, although the 
mechanism for such a deeper temperature rise is unknown. 
Alternatively, the discrepancy could be because the scatter- 
ing species is higher in the atmosphere than our MgSiOs 
example case. Our pressures can shift by an order of magni- 
tude assuming an MgSi03 haze, and could shift even more if 
another species is responsible for the scattering continuum. 

It is important to note that, if our pressure scale shifts 
significantly, our assumption that H2 is the dominant at- 
mospheric species may become invalid and, in this case, the 
2800 K temperature would be an overestimate. Conversely, if 
the abundance decreases significantly over the ranges where 
we measure the individual temperatures, for example due to 
ionisation, this will cause us to underestimate the temper- 
ature. Further theoretical work is required to better under- 
stand the overall properties of the atmosphere, especially as 
the models for the upper atmosphere of HD 209458b may 
not adequately describe HD 189733b. 



4.2.2 Comparing the T-P Profiles of HD 209458b and 
HD 189733b 

Fi gure[T2lalso shows the T-P p rofile derived for HD 209458b 
bv iVidal-Madiar et al.l <|2011ah , which was calculated using 
the same method as in this paper, of fitting temperatures 



to the spectral absorption depth profile . It is an update to 
the pre vious profile for HD 209458b bv iVidal-Madiar et al] 
(|2011bl ). which was determined by fitting temperatures to 
the integrated absorption depth profile, where measure- 
ments of the slope of the line can be compromised by di- 
lution effects of increasing bandwidth around an unresolved 
line. The authors find that the two analyses are equivalent 
for regions where the line is resolved and firmly detected, 
and their conclusions remain unchanged, but the error anal- 
ysis is non trivial due to all of the points in such a pro- 
file being co-depen dent. For comparison t o the HD 209 458b 
data, the mode l s oflShowman et all (|2009l 1. lYelle] <|2004h and 
iGarcfa Mufiozl (I2OO7F " are shown. Comparing the calculated 
T-P profiles for both planets to the models shows that the 
model profiles predict pressures lower than those derived 
from the observations even for HD 209458b at the ther- 
mobase. There are also no equivalent models for the up- 
per atmosphere for HD 189733b, since existing models use 
solar-type stellar spectra. The models of lMoses et al.1 (|201lh 
make allowances for the reduced emission at wavelengths 
longer than 2830 A, but the thermosp here is placed a t a spe - 
cific pressure, based on the models of I Gar cf a Mufiozl (|2007f ). 
Since there is no work at present that is able to model all 
of the factors at the base of the thermosphere or at low 
pressures for HD 189733b, we cannot draw any meaningful 
comparisons with existing models. However, it is hoped that 
these observations will be able to constrain future theoreti- 
cal work. 

The most useful comparison to draw currently is 
between the two T-P profiles calculated from obser- 
vations. To compare the two planets, we assume that 
our observed temperature rise indicates a detection of 
the base of the thermosphere, although we cannot 
be sure that the temperature continues to rise above 
our observed altitude regions. Both planets have mea- 
sured escaping atmospheres ( Vidal-Madiar et ail 1 20031 ; 



ILecavelier des Etangs et alj I201C ) , which should be associ- 
ated with extremely high temperatures at high altitudes. 

The shape of the derived T-P profile for the upper most 
regions of HD 189733b looks similar to the HD 209458b pro- 
file. However, this conclusion has to be tentative as the in- 
frared data points are not tied to our pressure scale, so the 
overall T-P profile over all observed pressures would look 
different if the upper atmospheric region of the profile were 
shifted. It would be very interesting to resolve the spectral 
absorption depth profile further and determine more accu- 
rately the shape of the T-P profile. The speed with which the 
upper atmosphere heats up with increasing altitude could 
give clues about any absorbing or reflecting material high in 
the atmosphere. 

4.3 The Relative Atomic Sodium Abundance 

Abundance measurements of identified spectral features can 
be determined from transmission spectra by rearranging 
Equation @. The relative abundance of two species is given 
by (also shown bv lDesert et al.ll2009T l: 



(6) 



§1 - £2 {(*i-*2)/*0 

where £1 and £2 are the abundances, o\ and 02 are the wave- 
length dependent cross-sections at Ai and A2, and Z\ and 22 
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Figure 12. Left: T-P profile for HD 189733b based on the assumption that the Rayleigh scattering signature is due to high-altitude 
atmospheric haze at a reference pressure of 10 — 4 bar. The temperatures are plotted in black for each pressure range (thick vertical 
bars). The thick vertical bars do no t include the uncertainty in refere nce pressure. The data include the 1340 K temperature and 



pressure derived for the continuum by Lccavclicr des Etangs et al. combined with the 1280 K temperature measured for the near 

infrared continuum by andlSing et al j 1 20091). The pl ot also shows lower altitude poin ts fromlKnu tson et al,l ||2007|) ; iKnutson et a.1.1 d2009fl . 



ICharbonneau et al.l d2008l ) and Deming et alj J2006T) ■ The temperat ures measured bv lCharbonneau et alj ll2008ri and lDeming et alH2006h 
were adjuste d to the terminator t emperature bv lHeng et al.l ll201ll) . T-P model profiles are shown for the lower altitude regions with red 
dashed lines jFortnev et al.ll2010h and for the high er al t itudes we show models for the upper a tmosphere of HD 209 458b for comparison, 
where the black solid line is the T-P profile fro m I Ye lle (2004) and t h e dot-d ashed line is from Garcia Munoz (2007). Right: T-P profile 
for HD 209458b using the measur ements fro m Vidal-Madi ar et al. I feOllah . plotted as black points. Show n with red dashed lines are 
mod els for the lowe r atmo sphere l|Showman et alj|2009l ). The solid black line is a model from lYelld (20041) and the dot-dashed line is 
from lGarci'a Munoj (120071) . 



are the altitudes of the observed spectral features at Ai and 
A2 of species 1 and 2 respectively. All other constants cancel 
out assuming that both spectral features occur in the same 
temperature regime. 

We compared our sodium absorption depth with the 
radius observed at 8 urn by lAgol et al.l (|2010l ). assum- 
ing that this feature is due to absorption from atmo- 
spheric water. The results from secondary eclipse observa- 
tions h ave detected a water absorption signal at this wave- 
length (|Grillmair et al. 2008), leading us to believe that wa- 
ter is likely to be the source of the 8 um feature w e see i n 
transmission as already envisaged by iDesert et all (2009). 
We used Equation (0 along with cth,o = 2 x 10 -20 cm 2 



at a wavelength of 8 urn ( Desert et al.l 120091 ) and a ra- 
dius for the wat er feature at 8 urn of R^/R* = 0.15531 
l|Agol et alj|2010l ). This was compared to the radius in the 
region 10-18 A from the centre of the Na I doublet, which 
was R-p/R* = 0.1572 ± 0.0002. This wavelength region has 
a similar temperature to the continuum and has absorp- 
tion cross-section oxa = 3.5 x 10~ 20 cm 2 . Since this is not 
a differential measurement, we used only the fitted radius 
from visit 7, which is has the smallest contamination from 
occulted stellar spots of all the visits, rather than an av- 
erage of all 3 visits. The un-occulted spot correction used 
for the G750M band was 1 per cent, and the un-occulted 
spot correction u sed for the 8 um band was 0.2 per cent (see 
ISing et al.|[201l| ). 

The solar abundance ratio of sodium to water is 
W£h 2 o Mar] = 6.63-9. 9 5X 10^ or ln(W£H,o [so l ar]) = 
-5.0 to -4.6 (iLoddersl 120031 ; lLodders fc Feglevl |2002| ; 



ISharp fc Burrows 
values of £Na/£H 2 o[189] 



20071 ). For HD 189733b, we calculated 
which were ~ 100 times higher 



(ln(^Na/^H 2 o[189]) = +1.3 ± 0.75), indicating a super-solar 
abundance ratio. 



5 SUMMARY AND CONCLUSIONS 

We have detected sodium in the atmosphere of HD 189733b 
with 9a confidence, confirming the previous ground-based 
measurement, and have improved the precision of the mea- 
sured absorption level by almost a factor of 3. Our improved 
resolution has enabled us to measure the sodium doublet ab- 
sorption profile. We confirm the presence of a very narrow 
Na I doublet feature, which can be explained by a high- 
altitude haze obscuring broad Na, K and H2O features, or 
a low Na I abundance hiding the wings beneath an H2 scat- 
tering signature. 

We use the Na I spectral absorption depth profile to de- 
termine the temperature as a function of altitude and find 
that two temperature regimes are required to explain the 
observed profile. The temperature rises with increasing alti- 
tude over the regions we have measured, indicating a likely 
detection of the planet's thermosphere. A thermosphere is 
expected from models and has been detected in HD 209458b. 
So far, only two planets have measured upper-atmosphere 
temperature profiles, but it is reasonable to assume that 
most hot Jupiters will exhibit hot thermospheres, due to 
their closeness to their host stars and the amount of radia- 
tion they receive in their upper atmospheres. 

We calculated a resulting T-P profile assuming a high- 
altitude haze as the reference level, based on MgSiOs as the 
scattering species. It is possible that condensates such as 
MgSiOs could sublime at temperatures of ~ 1300 K or even 
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lower tempera tures at the low pr essures sensed with the vis- 
ible data (e.g. iMoses et al1l201lf ). which could indicate that 
another species is responsible for the observed signature. 
The pressures are 3 orders of magnitude higher assuming H2 
Rayleigh scattering as the reference level, indicating how un- 
certain the pressure scale is. Additionally, the pressure scale 
with altitude will change if the sodium abundance is not 
solar. 

We find that we cannot constrain the Na I abundance 
without knowing the species responsible for the continuum 
signature and hence our reference level, but comparison with 
a feature at 8 urn indicates that the relative abundance of 
sodium compared to water is much greater than solar. 

We have measured only the second upper-atmospheric 
temperature-altitude profile for an extrasolar planet. The 
method described in this paper could be applied to other 
planets, leading to further understanding of their upper at- 
mospheres and how they interact with their host stars. The 
line profile gives information about whether high-altitude 
optical absorbers are present, allows us to determine an 
atmospheric temperature structure and will be invaluable 
for studying the atmospheric escape mechanism. It can 
also allow us to place some constraints on the elemental 
abundances if a reference pressure can be determined. It is 
thought that there are at least two classes of hot Jupiters, 
but more measurements are required to determine whether 
HD 189733b and HD 209458b are representative of the dif- 
ferent classes. The case of WASP-17b, w hich has an absorp - 
tion profile that is yet different again ()Wood et al.ll20l"lf l. 
suggests a range of properties within these classes, with an 
even narrower line profile suggesting an optical haze appear- 
ing higher in the atmosphere than in HD 189733b or an 
atomic sodium abundance lower than in HD 189733b. With 
regard to HD 189733b specifically, an optical secondary 
eclipse measurement would show us whether there is reflec- 
tion or scattering from absorbers high in the atmosphere, 
which will help constrain the nature of such absorbers. Ad- 
ditionally, high resolution ground based spectra are required 
to resolve this inner-most line core region, and extend the 
T-P profile to lower pressures and higher altitudes. 
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